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logical approach, we probe vacuum oscillations plus quantum decoher- 
ence between two neutrino species in the channels v„ — ► v T , v^ — » v e 
and z/ e — » v T , admitting that the quantum decoherence parameter 7 is 

qj ' related to the neutrino energy E v as : 7 = 70 E™, with n = —1, 0, 1 and 

2. Our bounds are valid for a neutrino mass squared difference com- 
patible with the atmospheric, the solar and, in many cases, the LSND 
K> . scale. We also qualitatively discuss the perspectives of the future long 

H I baseline neutrino experiments to further probe quantum dissipation. 
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I. INTRODUCTION 

The striking results of solar and atmospheric || neutrino experiments testify, 
beyond any reasonable doubt, that neutrino physics involves quantum interference 
phenomena. This is why it is plausible to envisage today the use of neutrino oscil- 
lations to probe the foundations of quantum mechanics and, in particularly, to test 
the completeness of the theory. 

From the theoretical point of view, quantum dissipative effects can be viewed 
either as a consequence of a fundamental violation of quantum physics 0,f§], moti- 
vated by quantum gravity, or of an effective description of a micro system weakly 
coupled to the macroscopic world in a open quantum system framework 0]. 

In the former approach, it is argued that quantum fluctuations of the gravita- 
tional field can lead to a loss of quantum coherence, making time-evolution transform 
pure quantum states into mixed ones, thereby violating ordinary quantum mechan- 
ics 0-§]. In the latter approach, quantum mechanics is not violated at the level of 
the global system but rather only by the reduced effective dynamics of the micro 
subsystem weakly coupled to the macro reservoir [H|§. 

Since both of these attitudes, although conceptually very different, will result 
in a modification of the time evolution of the neutrino mass eigenstates leading to 
the appearance of damping factors in the neutrino oscillation probabilities, we will 
not advocate in favor of either physical interpretation but rather simply treat here 
the effect phenomenologically. We consider that the absence of a full dynamical 
theory that can account for the origin, define the energy dependence and ultimately 
estimate the size of the decoherence effect is a good motivation for phenomenological 
analyses which can derive, directly from experimental data, limits which are valid 
independently of the (unknown) theoretical picture. 

Bounds on dissipative parameters have already been derived from studies of neu- 
tral meson systems PJlO Hllfl and neutron interferometry JL4]. A first attempt to use 



neutrino systems to investigate quantum dissipative effects was reported in Ref. IIS] 



In Ref. |16| , the possibilities of future neutrino experiments to unravel quantum de- 
coherence was qualitatively discussed. Recently, tight limits were obtained in the 
channel v^ —>■ v T from the atmospheric neutrino data, for a decoherence parameter 
7 which is either supposed to be a constant or proportional to the neutrino energy 
squared, and in fact a solution to the atmospheric neutrino problem (ANP) was 
found, when 7 was suppose to be inversely proportional to the neutrino energy |17| . 
The quantum decoherence parameter at the best-fit value of this novel solution to 
the ANP is such that it is not in conflict with other experimental constraints, but 
it is big enough to explain all the atmospheric neutrino data (sub-GeV, multi-GeV 
and upward going muons) comparably well as the mass induced v^ — > u T oscillation 
mechanism in vacuum. 

In this paper, we investigate the most stringent constraints on the quantum dis- 
sipation parameters one can get from terrestrial neutrino experiments, considering 
flavor conversions between only two neutrino species and assuming that these pa- 



rameters have some specific energy dependence. We will suppose that, in nature, 
mass induced vacuum oscillations are accompanied by quantum decoherence. In this 
framework, we will extract our limits on decoherence using data from experiments 
that have not registered any signal of flavor conversion, so that our bounds will be 
valid from a maximal value of Am 2 , which will depend on each experimental setup, 
down to Am 2 = 0. Although our bounds are more general, they will be true, in par- 
ticular, in the range of Am 2 consistent with the v^ — > v T oscillation solution to the 
ANP, i.e. Am 2 ~ (2 — 5) x 1CT 3 eV 2 [Tj| and the v e — > z/ M , v T solutions to the solar 
neutrino problem (SNP) either in vacuum, with Am 2 ~ 10~ 8 eV 2 — 10 -11 eV 2 |19|| , 
or in matter through the MSW mechanisms, with Am 2 ~ (2 — 20) x 1CT 5 eV 2 (LMA) 
or Am 2 ~ (4-10) x 10" 6 eV 2 (SMA) or Am 2 ~ (6-20) x 10" 8 eV 2 (LOW) 0. In 
the majority of cases for v^ — > v e , our constraints can also be applied in the LSND 
allowed region, 3 x 10~ 3 < sin 2 26 < 3 x 10" 2 and 0.2 eV 2 < Am 2 < 2 eV 2 [0-||. 



We will stress our results in these Am 2 regions, since they are, from the point of 
view of evidences in favor of flavor change, the most interesting ones. Moreover, we 
have to emphasize that in the point sin 2 Id = 1, our limits can be re-interpreted 
as a bound on flavor conversion driven by the pure decoherence mechanism alone 
(Am 2 = 0). 

After establishing the best current constraints, we briefly discuss the capability 
of future long baseline neutrino experiments to expand these limits. 

The outline of the paper is as follows. In Sec. |J, we review how to introduce 
quantum decoherence in the evolution equation of the neutrino mass eigenstates 
in the density matrix formalism. Using this modified formalism, we introduce the 
quantum decoherence parameters and justify the form of the neutrino oscillation 
probability we will use in this work. In Sec. |T|, we describe how we have analyzed the 
experimental data from CCFR |M|, E776 g§, CHORUS |^3TJ and CHOOZ jHJ 



in order to extract our limits. In Sec. [TV], we discuss the most important constraints 
from terrestrial experiments on the quantum decoherence parameters in the neutrino 
oscillation channels v^ — > u T , v^ —>■ v e and v e —>■ v T . In Sec. [V], we argue on the 
perspectives of KamLAND and other future neutrino experiments to put even more 
restrictive bounds on quantum dissipation. Finally, in Sec. 0, we present our 
conclusions. 



II. REVIEW OF THE FORMALISM 

The time evolution of neutrinos created in a given flavor v a by weak interactions, 
as of any quantum state, can be described using the density matrix formalism by the 
Liouville equation. In this formalism, the neutrino state in time can be described 
by a density matrix p a , which is a hermitian operator, with positive eigenvalues and 
constant trace. One can suppose that p a also has the additional property of being 
completely positive according to Refs. P,|32|, but this important theoretical point 



will not be crucial in the limit we are interested here. We will comment more on 
this below. 



Considering two neutrino generations, in the basis of the two mass eigenstates 
v\ and v-i that have masses ra\ and m 2 , respectively, the two flavor eigenstates v a 
and vp can be represented by 2 x 2 matrices as 

/ cos 2 9 cos 9 sin 9 \ , , 

Pa ~\cos9sm9 sin 2 # ' ^ 



and 
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sin 2 # 
— cos sin 9 



cos # sin 9 
cos 2 9 



Po 



(2) 



where 9 is the usual Cabibbo-like mixing angle that parametrizes the matrix which 
relates mass and flavor eigenstates. 

If we add an extra term L[p a ] to the Liouville equation, quantum states can 
develop dissipation and irreversibility |6|||. The generalized Liouville equation for 
p a {t) can then be written as || 



dtp a (t) = -t[H, Pa (t)] + L[ Pa (t)], 

where the effective hamiltonian H is, in vacuum, given by 
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(3) 
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where A = (mf. — m 2 )/4E„, we have already considered ultra-relativistic neutrinos 
of energy E„ and the irrelevant global phase has been subtracted out. 

One can rewrite Eq. (H) using as basis the identity and the Pauli matrices (a M , 
/i = 0, 1, 2, 3), here we drop the neutrino flavor index for simplicity, 



dtP^i Vfi = Ztnush^PvCTs + L^ u p v a^, 
where sum over repeated indices is implied, H = h^a^, p = p^a^ and 
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is the most general parametrization for L[p], it contains six real parameters which 
are not independent if one assumes the complete positivity condition |32] . While the 
requirement of simple positivity guarantees that the eigenvalues of p, interpreted as 
probabilities in the formalism, remain positive at any time, the requirement of com- 
plete positivity guarantees the same thing for the density matrix which describes the 
larger system (neutrino plus environment). This assures the absence of unphysical 
effects, such as negative probabilities, when dealing with correlated systems PJ3T 



The relations that must satisfy these parameters, if one assumes complete positivity, 
can be found in Ref. ||. 

Using the above parametrization, Eq. @ can be written as 

d t po{t) = 0, (7a) 

d t p x {t) = 2e ilvl h Vi p v - 2[opi + bp 2 + cp 3 ], (7b) 

d t p 2 (t) = 2t iiv2 h ii p v - 2[6pi + 7p 2 + PP3], (7c) 

dtP3(t) = 2e iiv3 h ii p v - 2[cpi + (3p 2 + ap 3 ]. (7d) 

These equations lead, in general, to a v a — > v$ transition probability that will 
present a time behavior characterized by two kinds of contributions: an oscillating 
one, controlled by A, and an exponentially damping one, driven by the dissipative 
parameters. The importance of these effects will depend on the relative magnitude 
of A with respect to the decoherence parameters. In particular, an asymptotically 
oscillatory behavior will be possible if these decoherence parameters are all small 
when compared to A. One can also envisage, as an opposite extreme case, the 
situation where decoherence is the dominant phenomenon. In this case there will be 
no oscillation but only a damping effect. 

Here we are interested in investigating the possibility of extracting limits on 
decoherence admitting that neutrino oscillation is the dominant process. The 
simplest situation is the one in which we can have the oscillation probability in 
the presence of a single extra dissipation factor. In fact, this situation physi- 
cally arises when the weak coupling limit condition is satisfied, and corresponds 
to a = — > a = 7, b = c = [3 = 0, if we impose complete positivity. The weak 
coupling limit is one way of implementing the physical requirement that the inter- 
action between neutrinos and environment is weak, as explained in Ref. || . Even if 
complete positivity is not assume, it is reasonable to think that all these parameters 
should take very small values, otherwise decoherence would be a predominant effect, 
so that the simplest situation would still be to consider that only one of them is 
large enough to be within experimental reach. In this case, Eqs. (0) simplify to 

Po{t) = p (0), (8a) 

dtfh (t) = 2 [h 2 p 3 - h 3 p 2 ] - 2 7Pl , (8b) 

d t p 2 (t) = 2[h 3 p 1 - hip 3 ] - 27p 2 , (8c) 

d tP3 (t) = 2[h lP2 - h 2Pl ]. (8d) 

Solving Eqs. (H), with hi — h 2 — and h 3 = —A from Eq. @ and the initial 
conditions which correspond to p Q (0) = p a , 

p (t) = i, (9a) 

Pi(t)= isin2#e- 2 ^cos(2A£), (9b) 

p 2 (t) = - | sin 20 e- 2 ^ sin(2At), (9c) 

p 3 (t)= i cos 20, (9d) 



so that the flavor conversion probability, which in any case can be calculated as 

P(u a ^u $ )(t)=Tr\p a (t)p $ (p)] 

= 1 - [i + p 3 (t)] cos 2 9 - [~ - p 3 (t)] sin 2 - Pl (t) sin 20, (10) 

gives in the simplest, but physically meaningful, limit we are interested in 

P{y a -> up) = - sin 2 29 [1 - e~ 27L cos(2AL)], (11) 

the probability of finding the neutrino produced in the flavor state v a in the fla- 
vor state vp after traveling a distance L under the influence of quantum dissipation 
driven by the parameter 7. If 7 = 0, we recover the usual mass induced oscillation 
probability in two generations. Although other possibilities exist, where all deco- 
herence parameters survive, this simple form will at least allow us to obtain from 
experimental data a very good idea of the order of magnitude of limits we can get 
on decoherence parameters. 

Notice that, by Eq. (|10|), even if there is no mixture, i.e. cos 9 = 1 and sin# = 0, 
there can be transition as long as ps(t) is not a constant, so that even if neutri- 
nos are degenerate or massless, flavor change can still take place through the pure 
decoherence mechanism (PDM). The simplest case of flavor conversion via PDM 
(Am 2 = 0) is the one where only a single decoherence parameter that appears in 
Eq. (7d), say a = 7, is non zero. This implies the probability of conversion 

P(v a -+v ) = \[l-e-*' L ]. (12) 



III. DESCRIPTION OF THE EXPERIMENTAL DATA ANALYSES 

We analyze data from the neutrino accelerator experiments CHORUS, CCFR, 
and E776, and the reactor experiment CHOOZ; which can currently put the most 
stringent limits on the dissipation parameter 7, in the context of two neutrino gen- 
erations. The exact dependence of the dissipation parameter 7 with the neutrino 
energy is unknown. Several forms, based on physical considerations, have been pro- 
posed in the literature ||[9 , [16 , 17 1 . For this reason, we will investigate in this work 



four different ansatz on the 7 energy dependence: 

7 = 70^, n= -1,0, land 2, (13) 

where E v is the neutrino energy in GeV and 70 is a constant given in units of 
GeV^ 1 -"). 

In the case of the specific baseline and energy range of the investigated experi- 
ments, with the mass squared differences that we will consider here, the cosine term 
in Eq. ( O) turns out to be 1, so that one can simply write 



P(u a ^v ) = -sin 2 29 [I -e^ 1 }, (14) 

where 7 is given by Eq. (|13D and (a,/3) = (fi,r), (//, e), (e, r). The survival proba- 
bility P[y a — > z/ a ) = 1 — P(y a — > z/g), in any given channel, merely by conservation 
of probability Note that, for sin 2 29 = 1, the above formula is equal to the one for 
the PDM given in Eq. (|12|) , so that at this particular point, our results can also be 
interpreted as a direct limit on pure decoherence. In spite of that, the two situations 
are not physically equivalent. 

From Eq. (|J"4j), one can have a quick idea about the bound on 70 that can be 
reached by a given experiment. Due to the fact that we are dealing with experiments 
that give no evidence of neutrino conversion, the condition P{y a — > vp) < P min must 
be satisfied, where P m j n is the minimum probability (or sensitivity) that could be 
measured at each experiment. This means that 

[l--]<|Pi, («> 



thus 



In (l - 2 



P 



min 



„, <r V sin 2 29 J .. in -i9 n^K 1 -™) na\ 

70 ~ " (L/km) < E u /GeV >" X 10 GeV ' (16) 

where n = —1,0, 1 and 2, L is the neutrino flight length in km and < E u > the 
neutrino mean energy in GeV. This equation permits us to know in advance which 
experiment will put the best bound on 70 for a given n, in each oscillation mode. 

We derive constraints on 70 from three oscillation channels, v^ — > v T , v^ — > v e 
and v e — > v T : 

1. For v^ — > u T1 we focus our attention in the range of sin 2 29 consistent with 
ANP, i.e. 0.8 < sin 2 29 < 1.0. In all the cases of n, we work with CHORUS, 
our constraints being valid for Am 2 < 0.4 eV 2 . 

2. For v^ — > z/ e , we choose to stress two different regions of sin 2 29: 0.5 < sin 2 29 < 
1.0 and 6 x 10~ 3 < sin 2 29 < 6 x 10~ 2 . The first, we will call the large amplitude 
case (LA) and the second the small amplitude case (SA). Our interest in these 
regions is because they cover a similar range in sin 2 29 as the standard oscil- 
lation solutions to the SNP []. We remark that, LA encompasses the vacuum, 
the MSW large and the MSW low mixing angle solutions, while SA includes 
a big part of the MSW small mixing angle solution to the SNP and of the 



region allowed by LSND ED, KARMEN |22| and Bugey El experiments. For 



1 Note that since we are working in two generations Piy^ — > v e ) — P(v e -^ v ti) by T 
invar iance. 



n = — 1, we obtain the best constraints on 7 with CHOOZ for LA, and E776 
for SA, valid from Am 2 < 1 x 1CT 3 eV 2 and Am 2 < 5 x 10~ 2 eV 2 , respectively. 
For n = 0,1,2, CCFR gives the best limits in both cases, LA and SA, which 
are all valid from Am 2 < 2 eV 2 . 

3. For u e — > z/ r , we emphasize only the LA case, based on the same motivations 
explained for the v^ — > v e mode. We get the best bounds on 70 for n — — 1 with 
CHOOZ (Am 2 < 1 x KT 3 eV 2 ), for n = 0, 1 with CHORUS (Am 2 < 4 eV 2 ), 
and for n = 2 with CCFR (Am 2 < 20 eV 2 ). 

At the following sections, we will describe the data analysis performed for each 
experiment in order to calculate limits on 70 with a precise statistical significance. 

A. CCFR 



The CCFR experiment at Fermilab [p4] is a neutrino oscillation appearance 



experiment which provides stringent limits on three different oscillation channels: 
(1) v^ — > u T , (2) z/^ — > v e and (3) v e — > v T . The neutrino beam in this facility 
consists of about 98 % (z/ M + z/ M ) and 2 % [y e + z/ e ). The experiment has a 1.4 km 
oscillation length of which 0.5 km is the decay region. 

In order to extract the CCFR bound, we have performed a statistical treatment 
of the data similar to the one prescribed in Ref . |33j , that is, we have minimized the 
X 2 function 

15 / a robs A/"theo\2 

X 2 = E— T^^i (17) 

i=l °i 

where N° hs and <Tj are the experimental points and errors for the i-th bin, respec- 



tively, read of from the electron spectrum given in Ref. p4| . This spectrum has 15 
energy bins that go from 30 GeV to 600 GeV in visible energy and we assume here 
this to be equal to the neutrino energy. JV 4 theo is the theoretical prediction of the 
i-th bin under a certain oscillation hypothesis. Explicitly, 

N? eo U = (N cxp {(P(p e -► v e ) + 5 a3 P(u e -. v T ) m(E)) 

+ R(E)[6 a2 P(^ -> v e ) + 5 al P{v^ -+ v T ) m(E)}}}, (18) 

where iV exp is the expected flux of v e calculated without oscillations taken from Ref. 
p|j , R(E) is the ratio between the fluxes of v^ and v e [21|, m(E) is the probability 



of misidentifying v T as v e (18 %) times the ratio of v T to v e charged current cross 
sections |53| and a = 1, 2, 3 mean respectively channels (1), (2) and (3) as explained 
above. Also, () means that we have averaged the neutrino propagation length over 
the neutrino decay pipe and the neutrino energy over the bin width. We were able to 
reproduce the experimental exclusion curves for each channel reasonably well. We 



show our results for the three oscillation modes inspected by CCFR in Fig. [T] (a), 
they can be compared to Fig. 5 of the last paper in Ref. [Q . 

To extract the quantum decoherence limit for each oscillation mode, we use the 
following method : at fixed sin 2 29 we minimize the \ 2 given in Eq. ( |TTD as a function 
of 70, thus obtaining the bounds using \ 2 ~ Xmin = 2.70, 3.84 and 6.63 at 90 %, 
95 % and 99 % C.L., respectively. We have verified for each sin 2 29 that Xmin occurs 
at 7o = 0. 



B. E776 



The E776 [25] experiment at the Brookhaven National Laboratory searched for 



v ii, — > v e oscillations, looking for v e appearance 1 km from the source of a wide-band 
v^ beam of average energy around 1 GeV. 

The two primary backgrounds for this experiment were v e contamination in the 
beam and ^-induced 7r°'s which are misidentified as electrons by the detector. The 
E776 data collected was consistent with the expected background, so no signal of 
V/i -^ v e oscillation was found |25| . 

To reproduce the E776 limits on v e (z/ e ) appearance, we have minimized the sum 
Xtot = X 2 +X 2 ) which can be derived from a likelihood function and takes into account 
both Poisson statistical distribution of data and the presence of background |35,36|j, 
using the definitions 
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x 2 - 
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iV h + Br 



(19) 



where iV° bs [N i ) and A^ th (iV^ ) are the observed and the theoretical contents of 
v e (v e ) events in the z-th bin, in the presence of the background level Bi (Bi). 

N° hs (N° bs ) and B { (Bi) were taken from Fig. 3 of Ref. 
calculated as 



■t.b , 



2|] and A^f (N { ) 
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dE, 



(20) 



where the integration was performed over the z-th bin of the z/ M (i/^) spectrum N v 
(N v ), also given in Ref. p5| . The cross sections were taken from Ref. [p4 ]. We show 



our exclusion plot for the two generation mass induced oscillation in Fig. [I] (b), it 



agrees well with the E776 plot in Ref. |25 



To obtain the decoherence limits, we have minimized Xtot> an< ^ computed the 
excluded regions for 70 at various confidence levels, in the same way as already 
described for CCFR. 



C. CHORUS 



CHORUS p6|-p9| experiment at CERN is a neutrino oscillation experiment that 
has been looking for u T appearance in the channels v^ — > u T and v e — > v T . The 
CHORUS detector is illuminated by a beam consisting mainly of v^ coming from 
the "West Area Neutrino Facility". The average z/ M momentum is 27 GeV and the 
average distance they have to travel from the target to CHORUS is about 0.6 km. 
The u^ beam is also composed of 5% v^ and about 1% v e + v e . The genuine v T 
content of the beam is estimated to be negligible, at the level of 3.3 x 1(U 6 v T 
charged current interactions per v^ charged current interactions ||28|| , well below the 



sensitivity limit of the experiment. The aim is to detect v T by observing the r~ 
produced in a charged current interaction and its subsequent decay vertex in an 
active nuclear emulsion target. 

The basic difference between the z/ M — > v T and v e —>■ v T searches comes from the 
energy spectra, the average energy of v e is about 20 GeV higher than u^, which leads 
to differences in the acceptances for v T interactions [ES . 



The main sources of background are the charm D meson production followed 
by its muonic decay (for the leptonic modes), the D meson hadronic decays and 
hadronic scattering on target nuclei without any visible nuclear breakup or recoil 
(for the hadronic modes) and prompt v T from the beam [^,^ . 

We define the average oscillation probability 

J$ Ua (E)P(v a ^v T )dE 

where a = /x, e, $„ and $„ e are the v^ and v e fluxes, respectively, taken from Ref. 
2(J and P[y a — > u T ) is the oscillation probability given in Eq. (|l"4]) . The exclusion 



region can be calculated according to Refs. pqj29|1 by using 

(Pa^r) < ^p, (22) 

where ^ xp involves ratios of cross sections, acceptances, efficiencies, number 
of muonic and non-muonic events and branching ratios for the different decay 
modes [ES| and F is a numerical factor for zero v T events observed, no background 



expected and a total systematic error of 17% following the prescription of Ref. |37] 
i.e. F = 2.38, 3.12 and 4.9 respectively for 90%, 95% and 99% C.L 



From Ref. p9| , we can extract: ^ xp = 5950.0 for ^ -+ v T and iV e exp = 79.3 for 



v e — > v T . Using these experimental inputs in Eq. (^) we were able, for 70 = 0, to 
reproduce quite well the exclusion regions at 90% C.L. presented by CHORUS. This 
can be seen by comparing our plot for CHORUS in Fig. |l|(c) with Figs. 1 and 2 of 
Ref. 0. 

The exclusion region for quantum decoherence was obtained by fixing sin 2 29 and 
computing the 70 value which satisfies each confidence level requirement given by 



Eq. (^). We did not do the analysis for the latest CHORUS result of Ref. |3(| since 

10 



we do not know exactly what is the total systematic error of this data, hence the 
correct value of F we should use to re-interpret their results. In any case, this will 
not change very much our conclusions. We will comment more on this latter on. 

D. CHOOZ 



The CHOOZ reactor experiment |yj is the one which provides the most re- 



strictive v e — > u e disappearance limit. The baseline is about 1 km and the neutrino 
energy around 3 MeV, so that L/E ~ 300 m/MeV. Since it is a disappearance exper- 
iment, it can be employed to infer bounds on 70 in the oscillation modes v e — > v^ and 
v e —>■ v T . In order to do this we have used the average of the ratio of the measured 
positron spectrum over the expected one, R = 1.01 ± 0.028(stat) ± 0.027(syst) [[H|] 
and compared this with the v e — > v e survival probability averaged over the spectrum. 
We work with the following x 2 definition 

^/q-pfr. -.,g)- ay (23) 

where a = |i 01 r, () means that the survival probability has been averaged in 
the positron energy spectrum and or = 0.039 is the total error in R. We recall 
that the positron energy and the neutrino energy are related by a simple rescale: 
E e + = E Ve — 1.804 MeV which allows us to easily perform the above mentioned 
average. The oscillation probability P[y e —>■ u a ) is given in Eq. (|14]) . Again, the 
excluded region is computed using the method already described in Sec. [ HI A| . We 



have adopted this elementary \ 2 formula since it allows us to reproduce almost 



exactly the exclusion curve obtained by the analysis A of Ref . |3Jj . This can be seen 
in Fig. [[] (d). We have checked that the global minimum of the y 2 always occurs at 
7o = 0. 



IV. CONSTRAINTS FROM CURRENT EXPERIMENTAL DATA 

Here, we present the limits on 70 given by experimental data in three different 
neutrino oscillation modes. We display in Figs. Q-^ our exclusion plots in the plane 
sin 2 29 x 7 at 99% C.L., they give a general view of our results and show the 
sensitivity of each experiment for each oscillation mode. 

Experiments with negative results on neutrino oscillations, such as the ones 
studied here, have a natural restriction to impose limits on 70. This happens when 
sin 2 29 — > 2P min . At this point the right hand side of Eq. (|T5| ) goes to 1, then for 
any value of 70 the left hand side (that can only go from (70 — > 0) to 1 (70 — > 00)) 
satisfies the inequality, so that we are not sensitive to variations of 70. Because of 
this fact, we cannot explore the SA case for u e — > z/ r , or use CHOOZ, instead of 
E776, to put restrictions on 70 in the mode v^ — > v e with n = — 1 in the SA case. 
These features can be seen in Figs. |3] and |[ 

11 



In the following, we will discuss the regions of special interest, in each mode, as 
underlined in Sec. fTT[ 



The best limits in this oscillation mode for all values of n are given by the 
CHORUS experiment, see Fig. ||. 

Previous to comment our bounds in this channel, we will make a brief remark 
about the constraints on 7 from atmospheric neutrino experiments, since they seem 
to imply Vp — > v T oscillations []38| . These experiments are clearly in great advantage 
over terrestrial ones, they cover oscillation baselines from approximately 500 km to 
about 12 000 km, the diameter of the Earth, with neutrino energies from ~ 1 GeV 
up to a few hundred GeV. One can estimate the order of magnitude of their bounds 
on 70 by asking 0{2^L) < 1 using L ~ 10 4 km and < E v >~ 10 2 GeV (the most 
energetic upward going neutrinos will push the limit), this gives : 70 < 10~ 21 GeV 2 
for n = -1, 70 < 10~ 23 GeV for n = 0, 70 < 10" 25 for n = 1 and 70 < 10~ 27 GeV" 1 
for n — 2. This is in good agreement with the limits obtained in Ref. [|I7| f\ for 



n = —1,0 and 2 using a statistically rigorous analysis of the atmospheric neutrino 
data. 

We display the bounds on 70 as a function of sin 2 29 in Fig. |2|, for n = —1, 0, 1 
and 2 at 99% C.L. In this figure the lower limit of the amplitude region compatible 
with the standard solution to the ANP is marked by a line with an arrow. For n = 2, 
the analysis of the CCFR data described in Sec. pi A] also can provide a very similar 
limit, slightly less restrictive than the one given by CHORUS. We do not show this 
here. 

We observe that some of our bounds are substantially weaker than the ones 
given by the atmospheric data. Notwithstanding, they are independent limits and 
the most restrictive ones that can be computed with data from neutrinos produced 
in accelerators up to this day. 

Finally, it is easy to re-scale our results to obtain more stringent constraints for 
the combined CHORUS/NOMAD actual data analysis |39|]. In fact, it is possible 



that in the near future CHORUS and NOMAD combined analysis may even provide 
a better bound than the atmospheric neutrino experiments for the case of n = —1. 
However, it seems unlikely that they will be able to overcome the atmospheric limit 
on 70 for n = 0, 1 and 2. 



2 We remark that in our notation 270 corresponds to 70 of Ref. |L7|] 
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B. Vfj, ->• v e 

The CHOOZ/E776 experiments provide the best limits for n = — 1 and CCFR 
for n — 0,1 and 2. This can be seen in Fig. |3|. Here, we focus on two extreme 
situations: LA and SA. 

Before discussing our bounds in this channel, we need to introduce a comment 
linking the SNP and quantum decoherence. This is in order to make a comparison 
between our results and the possible ones in the Sun. 

In the case of solar neutrinos, we can envisage two different approaches: (i) 
decoherence influences long wavelength oscillations or (ii) decoherence perturbes 
MSW-type oscillations. In the event of (i), L has to be taken as the Sun-Earth 
distance, L ~ 1.5 x 10 8 km, in the event of (ii) as the Sun radius, L & ~ 7 x 10 5 km. 
Using the fact that, the quantum decoherence parameter in order to be revealed 
must satisfy 0{2^L) ~ 1, and the range of neutrino energy is 10~ 4 GeV < E u < 
1CT 2 GeV, we can estimate the order of 7 for both situations. Considering (i), we get 
7o ~ 7x{(l(r 32 -10- 30 ) GeV 2 , 10" 28 GeV, lO^-lO" 24 , (l(T 24 -iO- 22 ) GeV" 1 } and 
considering (ii) 7o ~ 1.4 x {(10~ 29 - 1(T 27 ) GeV 2 , 1(T 25 GeV, 1(T 23 - 10~ 21 , (1(T 21 - 
10 -19 ) GeV -1 }, respectively for n = —1, 0, 1, 2. This means that the solar neutrino 
data can at the most set bounds on 7 within the ranges quoted above. 

In the case of LA, we show our results in Fig. |3] for n = —1, 0, 1 and 2 at 99% 
C.L., we see for n = —1, that the CHOOZ and CCFR limits are weaker than the 
possible constraints one can get from solar neutrinos. 

For n = l, the CCFR constraints are in the fringe of the solar sensitivity for 
long wavelength oscillations. Nevertheless, given that our results at sin 2 29 = 1 can 
also be understood as limits on the PDM alone, we can conclude that, for n = 1 
the possibility of explaining the solar neutrino data by this mechanism is in fact 
discarded. This is because not even the total rate measured by Super-Kamiokande 
can be explained if 70 < 3 x 10~ 24 , unless we admit that the 8 B flux can be 50 % 
smaller, i.e. 3 sigma away from the central value of the Solar Standard Model 
prediction. For MSW-type oscillations and n — 1, the CCFR limits are clearly more 
restrictive than those that could be reached by solar neutrino data. Obviously, the 
solar neutrino problem cannot be solved by the PDM in this case. 

When n = 2, the CCFR limits are much lower than the solar sensitivity for both 
(i) and (ii). Here again, the PDM cannot be solution of the solar neutrino problem. 

The general result is that for n > 2 (i) and n > 1 (ii), one cannot hope to 
improve the CCFR bounds with solar data nor to explain the SNP by the PDM. 
This is due to the energy dependence of the limit on 70 (see Eq. (|T3|)) and the fact 
that CCFR average neutrino energy is much higher than that of solar neutrinos. 
Therefore, the tendency is that the CCFR limits become stronger than the solar 
neutrino ones, with increasing n. 

For SA, we can also see the limiting curves in Fig. |3|. Note that, in this case, 
the constraints for n = — 1 were derived using E776 data. 

For the case n = —1,0 and 1, we expect that the bounds with solar neutrinos 
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should be much better than the ones we present here. In the case n = 2, the CCFR 
limits are below the solar sensitivity only for the situation (ii). 



As mentioned in Sec. |HI| , the SA bound for v^ — > v e mode can be kept also valid 
under the hypothesis that the oscillation parameters lay in the LSND/KARMEN/- 
Bugey HD^I allowed region (3 x 10" 3 < sin 2 2K3x KT 2 , 0.2 eV 2 < Am 2 < 
2 eV 2 ). Due to this range of Am 2 , just the bounds with CCFR data (n = 0, 1,2) 
can be used if the LSND mass scale is adopted. 

To compare our constraints with the possible ones that could be extracted from 
the LSND data, we have to do some estimations. Now, working with the assumption 
O^L) ~ 1 applied to LSND using < E v >= 40 MeV and L = 30 m, we arrive 
at 70 ~ {10~ 19 GeV 2 ,3 x 10~ 18 GeV,0.8 x 10~ 16 ,2 x lO-^GeV" 1 } respectively for 
n = —1, 0, 1, 2. This means that, for n — — 1, we are limiting the same range of 70 
that could be achieved with the LSND data. Meanwhile, for all the others n the 
CCFR limits are much below the LSND sensitivity. Furthermore, if we compare 
these sensitivities with our bounds on the PDM (sin 2 26 = 1), it is evident that the 
LSND results cannot be explained solely by decoherence. 



C. z/ e — » v T 

In this mode, the CHOOZ experiment will provide the best constraints when 
n = —1, CHORUS when n = 0, 1 and CCFR when n = 2, as can be seen in Fig. §. 

As we already mentioned in the beginning of this section, here we can only discuss 
the LA case. We verify in this oscillation mode that only for n > 2 the experimental 
limits imposed by CCFR are stiffer than those that could be achieved by the solar 
neutrino data. 



V. FUTURE PERSPECTIVES 

In this section, we briefly discuss the perspectives to improve the bounds 
presented before, using the mean values of E u , L and the sensitivities ex- 
pected for the forthcoming neutrino oscillation experiments such as Kam- 
LAND H, MINOS @-H, the CERN-to-Gran Sasso neutrino experiments ICA- 



NOE/OPERA [0,0 and a possible neutrino factory in a muon collider |46| . 

To estimate the limits on 70 that could be reached at prospective long baseline 
facilities, we investigate in turn the two possible outcomes of the experiments : (i) no 
signal and (ii) a positive signal of neutrino appearance/disappearance is observed. 

Experiments that give negative results are insensitive to the mass scale, so when 



considering (i) we can work with Eq. (I5|), replacing sin 2 28 by 1, and the mean 
value of E u and L as well as the sensitivity goals by the corresponding ones of each 
proposed experiment. In this manner, at this fixed value of sin 2 29, we can test 
simultaneously decoherence plus mass oscillation and the PDM alone. This is done 
for v^ — > v e and v e — > v T modes only, v^ — > v T is assumed to be already observed 
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by Super-Kamiokande and K2K. This will be a rough estimate since we completely 
neglect matter effects, which could be crucial in determining the exact sensitivity of 
each experiment to 70. Nevertheless, we believe that this will affect only the actual 
value of the limit but will not change our conclusions on which will be the most 
restricting experiment at each n and oscillation mode. 

In Fig. |5|, we plot our limit on 70 in the v^ — > v e and v e — > v T oscillation channels 
as a function of n for a few experimental configurations in case that no oscillation 
is observed. 

Let us discuss the v^ — > v e mode first. For n = — 1, KamLAND can bring the 
current bound three orders down, as well as certain configurations of a neutrino 
factory in a muon collider as exemplified in the plot. For n — 0, MINOS/ICANOE, 
since they have the same L and proposed sensitivity, can improve the limit by three 
orders of magnitude. OPERA will be less powerful since its proposed sensitivity 
is somewhat lower. A neutrino factory can bring this limit almost six orders of 
magnitude down. For n = 1, MINOS/OPERA/ICANOE can lower the bound by 
about a factor ten, but the neutrino factory can push it close to 70 < 10~ 29 . In the 
case n — 2, the CCFR constraint can only be overcome by a future neutrino factory. 

In the v e — > v T mode for n = — 1, KamLAND can certainly bring the CHOOZ 
limit down by a few orders of magnitude. On the other hand, some of the proposed 
neutrino factories can get better limits on 70, for all the values of n, than the ones 
discussed in this paper. 

When considering (ii) to put a bound on 70, the oscillation probability must be 
affected by the decoherence parameter so that 0(2 / yL) ~ 1. Thus, this value will 
not depend on the mode of oscillation, but only on the characteristic L and mean 
E v of each experiment. This gives us an idea of their sensitivity to 70 as a function 
of n which is shown in Fig. ^J for a few situations. 

For v^ — > v T only the ansatz n = — 1 can, in principle, be better tested in a 
neutrino factory than by using the atmospheric neutrino data. This will greatly 
depend on the specific choice of baseline and of the muon energy. 

In the mode v^ — > u e , the bound on 70 with the ansatz n = — 1 can be improved 
by KamLAND. For n > 2, the CCFR limits we give in this paper can be lowered 
by some of the experiments we studied. 

In the mode v e -► u T , KamLAND (n = -1), MINOS/ICANOE/OPERA/neutri- 
no factory (n = 0, 1) and a neutrino factory (n = 2) can improve the present limits, 
but only by one to two orders of magnitude. 

VI. CONCLUSIONS 

We have analyzed the experimental data from the terrestrial neutrino oscillation 
experiments CHOOZ, CHORUS, E776 and CCFR in order to extract constraints on 
the decoherence parameters in the framework of two generation neutrino oscillations 
induced by mass plus quantum dissipation. This was done for the so called weak 
limit, where all decoherence parameters, except for 7, are supposed to be either zero 
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or out of experimental reach. We have derived constraints on 70 under four different 
ansatz, i.e. 7 = 70 E™, with n = —1,0,1,2, from examining experimental data in 
the channels z/ M — > u T , v^ — > v e and v e — > v T . This is valid from distinct upper 
values of Am 2 , depending on each experiment so that we can always consider Am 2 
negligible in the probability, down to Am 2 = 0. In particular, we discuss our limits 
paying attention to the phase space of (Am 2 , sin 2 29) consistent with the standard 
solution to the ANP for v^ — > u T , and to the SNP in the channels v e — > v^ (LA, 
SA) and v e — > v T (LA). Our limits on 70 for v^ — > v e (SA) are also valid for Am 2 in 
the LSND allowed region. In addition, at sin 2 29 = 1, our bounds can be also read 
as direct limits on the PDM. At the end, we have also discussed the perspectives of 
future experiments to push our limits further down. 

In the Up — > v T mode (ANP solution range), we have established the following 
bounds: 70 < (5.6 - 4.3) x 10~ 21 GeV 2 for n = -1, 70 < (1.6 - 1.2) x 1(T 22 GeV 
for n = 0, 70 < (3.2 - 2.4) x 1(T 24 for n = 1 and 70 < (4.0 - 3.1) x 10~ 26 GeV^ 1 for 
n — 2, at 99% CL. In spite of the fact that these limits are much less restrictive than 
the ones given in Ref. [[17]] from analyzing atmospheric neutrinos, they are valuable 
to be known, since they are independent constraints. 

In the Up — > u e mode (LA), we have established the following bounds: 70 < 
(2.5 - 1.2) x 10~ 22 GeV 2 for n = -1, 70 < (6.0 - 3.1) x 10~ 22 GeV for n = 0, 
70 < (5.5 - 3.0) x 10~ 24 for n = 1 and 7 < (2.2 - 1.2) x lO" 26 GeV" 1 for n = 2, at 
99% C.L. From these constraints, we concluded that for n > 1 one is discouraged 
to try to extract better limits from the solar neutrino data itself. Moreover, these 
constraints exclude any possibility to explain the LSND results by the PDM alone. 

In the U/j, — > v e mode (SA), we have established the following limits: 70 < 
(6.0 - 0.27) x 10~ 19 GeV 2 for n = -1, 70 < (7.0 - 0.6) x 10~ 20 GeV for n = 0, 
70 < (7.0 - 0.5) x 10" 22 for n = 1 and 70 < (8.0 - 0.3) x 10~ 24 GeV" 1 for n = 2, at 
99% C.L. In the case n > 2, the solar neutrino data will give weaker bounds than 
ours. Besides, for n > 0, our constraints are stronger than what we could obtain 
with LSND data. 

In the v e — > v T mode (LA), we have established the following limits: 70 < 
(2.5 - 1.1) x 10~ 22 GeV 2 for n = -1, 70 < (1.0 - 0.5) x 10~ 20 GeV for n = 0, 
70 < (1.3 - 0.7) x 10~ 22 for n = 1 and 70 < (2.0 - 1.0) x 10" 24 GeV" 1 for n = 2 at 
99% C.L. Here again, for n > 2, our bounds are stronger than the solar ones. 

In the future, there will be many experiments capable of further constraining 
7q. However, it is rather difficult to say something completely definite about this 
since these experiments may or may not observe a positive signal of oscillation. 
In general, we can say that in the v^ —>■ v T mode only for the ansatz n = — 1 
the current limits from atmospheric data, can be improved by a future neutrino 
factory; in the v^ — > v e and v e —* v T channels KamLAND in = —1) and a future 
neutrino factory {n = 0, 1, 2) are the best candidates to test decoherence. We point 
out that in v^ — > v e mode (LA), for n > 2, the constraints discussed here can 
hardly be overcomed by any prospective experiment that observes a positive signal 
of oscillation. We remark that for n < 1 the solar neutrino data is certainly the best 
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probe for quantum decoherence in v^ — >■ v e and v e — > v T modes. 

Recently, bounds on pure quantum decoherence have been determined in Ref. Wi 



from SN1987A data. In fact, using our assumptions on the 7 energy behavior and 



the results of Ref. P8| |, we can estimate that SN1987A can provide the following 
limits, using L ~ 50 kpc and P(u e -> v^v T ) < 0.14 for (E„) = 26 MeV @: 
70 < 0.6 x 10" 39 GeV 2 for n = -1, 70 < 0.2 x 10~ 37 GeV for n = 0, 70 < 0.8 x 10~ 38 
for n = 1 and 70 < 0.3 x 10~ 34 GeV -1 for n = 2. It is very important to note that, 
these limits are not comparable to our bounds. This is because we are assuming 
decoherence plus mass induced oscillation with 10 -11 eV 2 < Am 2 < 2 eV 2 . In 
this context, the oscillation probability for large propagation distances will be no 
different than the one for mass induced oscillation, i.e. P{y a — > v$) — > sin 2 29/2 
(see Eq. (|i~T|) ), therefore independent of 7. Consequently, supernovae or other far 
away astrophysical objects emitting neutrinos cannot provide limits on the quantum 
decoherence parameter 70 in our framework. 

On the other hand, even in the situation that one can compare the SN1987A 
limits with the ones obtained here, i.e. for the PDM (sin 2 29 = 1), the constraints 
extracted from reactor and accelerators experiments, where the neutrino fluxes are 
well controlled, are very robust and worthwhile to be known. 
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FIG. 1. Our reproduction of the limits given by CCFR (a), E776 (b), CHORUS (c) 
and CHOOZ (d) on neutrino oscillations. The excluded region at 90% C.L. is the one to 
the right of each curve. 
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FIG. 2. Limits on 70 as a function of sin 29 in the v^ — > v T mode, for n = —1,0, 1 
and 2. The excluded region at 99 % C.L. is the one to the right of each curve. The lower 
value of the amplitude consistent with the ANP is marked by a dashed line with an arrow. 
All the limits were obtained with CHORUS data. 
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FIG. 3. Limits on 70 as a function of sin 2 29 in the i/„ — ► f e mode, for n = — 1 (SA, 
LA), 0,1 and 2 at 99 % C.L. 
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FIG. 4. Limits on 70 as a function of sin 20 in the u e 
and 2 at 99 % C.L. 
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FIG. 5. Limits on 70 for some of the next generation neutrino long baseline experiments 
assuming the proposed sensitivity of each experiment and no signal of oscillation observed 
in the v u — » v e and v e — ► v r modes. 
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FIG. 6. Attainable limits on 70 for some of the next generation neutrino long baseline 
experiments assuming they will observe neutrino oscillation. 
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